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ABSTRACT 

We report multiwavelength observations of the black hole transient GX 339-4 during its outburst 
decay in 2011 using the data from RXTE, Swift and SMARTS. Based on the X-ray spectral, temporal, 
and the optical/infrared (OIR) properties, the source evolved from the soft-intermediate to the hard 
state. Twelve days after the start of the transition towards the hard state, a rebrightening was 
observed simultaneously in the optical and the infrared bands. Spectral energy distributions (SED) 
were created from observations at the start, and close to the peak of the rebrightening. The excess 
OIR emission above the smooth exponential decay yields flat spectral slopes for these SEDs. Assuming 
that the excess is from a compact jet, we discuss the possible locations of the spectral break that mark 
the transition from optically thick to optically thin synchrotron components. Only during the rising 
part of the rebrightening, we detected fluctuations with the binary period of the system. We discuss 
a scenario that includes irradiation of the disk in the intermediate state, irradiation of the secondary 
star during OIR rise and jet emission dominating during the peak to explain the entire evolution of 
the OIR light curve. 

Subject headings: black hole physics — ISM: jets and outflows — X-rays: binaries 



1. INTRODUCTION 

Galactic black hole transients are binary systems that 
undergo sporadic outbursts that last for months. Dur- 
ing an outburst, they exhibit two main X-ray spectral 
states: the hard state and the soft state. In the soft 
state the X-ray energy spectra are dominated by a ther- 
mal disk component, and weak or no variability is ob- 
served; whereas in the hard state, the X-ray energy spec- 
tra are dominated by a non-thermal component, and high 
variability is observed. There also exists intermediate 
states in which the X-ray properties of the systems do 
not suit to the main states, bu t present the properties 
of a mixture of both states (see [Belloni 2 0101 for the de- 
tails of the spectral states). Besides the correlated X-ray 
spectral and temporal properties, black hole transients 
also show state dependent radio, optical and infrared 
(OIR) prope rties. In the soft state, the radio emis sion 
is quenched ([Fender et al.lll999t IRussell et aTll2011al) in- 
dicating the jet turning off. In the hard state, com- 
pact steady jets are observed (jFender! [2006( 1. Therefore, 
from observational point of view the transitions from the 
soft to the hard state serve the perfect conditions to 
study the properties of the accretion flow for jet forma- 
tion. Also in the hard state, the X-r ay flux is positively 
corre l ated with both th e radio flux dCorbel et al.ll2000L 
200 1 iGallo et~aT1[200l and the OIR flux ([Russell et al.l 
20061 ). These relations suggest a common origin between 
the jet emission and X-rays. Some jet models assume 
that a hot electron corona is the base of the jet to ex- 
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plain these relations (jMarkoff et al.l [20031 ). There are 
also reports of direct jet synchrotron emission dominat- 
ing the entire X-ray emission in t he hard state at 10 ~ 4 < 
^Edd < 10" 3 ([Russell et alJ2010t ). S oftening of the X-ray 
spectrum observed in some sourc es dTomsick et al.ll2001t 
Corbel et all 120061 : iDincer et all 120081: iSobolewska et~atl 



20111 ) at very low flux levels also suggests a change in the 



X-ray emission mechanism. 

The X-ray observations of GX 33 9-4 in January 
2010 r evealed the sta rt of an outburst ([Yamaoka et al.l 
l2010t iTomsickl 12010ft . Its multiwavelength observa- 
tions during the rise and spectral properties dur- 
ing the state transitions have been reported else- 
where (ICadolle Bel"eTall [20111 IGandhi et all l2?TTlt 



Shidatsu et all I2011al lbl: IStiele et all 12011b lYan fc Yd 
20111 ). In January 2011, the source made a transi- 
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tion from the int ermediate to the hard stat e during the 
outburst decay ([Munoz-Darias et al.l 1201 TT ) . The opti- 
cal observations throu gh the end of Februa ry 2011 re- 
vealed a rebrightening ([Russell fc lewisll2011[ ) which was 
als o observed in prev ious outburst decays of GX 339- 
4 ([Buxton et all 120121 ). Here, we report the results of 
RXTE, Swift, and SMARTS observations of GX 339-4 
in 2011 decay. We characterize the evolution of the X- 
ray spectral, temporal, and OIR photometric properties, 
and produce spectral energy distributions (SED) in order 
to investigate the jet-disk connection. 

2. OBSERVATIONS AND ANALYSIS 

2.1. RXTE Observations 

The outburst decay was amply covered with 54 pointed 
RXTE observations between MJD 55,560 and 55,650 
(2010 December 30 and 2011 March 30). However, some 
observations were not statistically satisfactory due to the 
short good time intervals (GTI) and/or small number of 
Proportional Counter Unit (PCU) on during the opera- 
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tion. Therefore, we did not include the observations with 
GTI < 500 s (see Table Q] for a log of observations) . 

We used data from the Proportional Counter Array 
(PCA) i nstrument onboard the RXTE for the spectral 
analysis (jJahoda et al.l [19961) . In most of the observa- 
tions, the spectra were extracted in the 3-25 keV energy 
band, but in a few cases for which the noise dominated 
above 20 keV, we used the 3-20 keV band. The response 
matrix, and the background model were created using the 
standard FTOOLS (v6.11) programs. We added 0.5% 
systematic error to the energy spectra following the sug- 
gestions of the RXTE team. 

The spectral anal ysis were performed using XSPEC 
12.0.7 ()Arnaudll 19961 ). We employed a spectral model for 
the continuum that consists of absorption, a multicolour 
disk black-body and a power law. We also included a 
phenomenological smeared edge model (Ebisawa ct al. 
[1991 for the iron K a absorption edge seen around 7.1 
keV to obtain acceptable x 2 values for the observations 
before MJD 55,606. In the spectral fits, the hydrogen col- 
umn density, Njj , and the smeared e dge width were fixed 
at 0.5 x 10 22 cm" 2 (jKong et al.ll2000l ) and 10 keV, respec- 
tively. This model was used previously in lTomsick et aLl 
(l2^0l : lKalemci et all (12004 12005L l2006aT) . 

The Galactic ridge emission was an important factor 
for the faint observations. In order to estimate its spec- 
trum we compared quasi-simultaneous RXTE and Chan- 
dra observations obtained on MJD 52,911. We combined 
seven RXTE/PCA observations taken on the same day 
and fitted with a model consisting of interstellar absorp- 
tion, a power law and a Gaussian to represent the Galac- 
tic ridge emission, and a second power law to represent 
the source. The centroid energy and width of Gaussian 
was fixed at 6.6 keV and 0.5 keV respectively. The pa- 
rameters of the second power law from the source was 
set to the yalues obtained from Chandra observation 
(|Gallo et al.ll2003D . With this method we modelled the 
Galactic ridge emission with a power law index of 2.1 
and an unabsorbed flux of 7.55 x 10 -12 ergs cm" 2 s _1 



in th e 3-25 keV band (jDincer et all 120081 ICoriat et all 
I2009j) . We applied Galactic ridge emission correction to 
the spectra, fluxes and rms amplitudes of the observa- 
tions after MJD 55,627. The contamination from Galac- 
tic ridge emission was not greater than 3% of the total 
flux for the observations before this date. 

For each PCA observation, we produced power den- 
sity spectrum (PDS) from segments of length 256 s with 
a time resolution of 2 -9 s in 3-30 keV energy band us- 
ing IDL programs developed at University of Tubingen. 
The averaged PDS were no rmalized as described in 
iMivamoto fc Kitamotol ( 19891) and co r rected for the dead 
time effects according to lZhang et al.l (|1995l ) . Then, PDS 
were fit using Lorentzians in 0.004-256 Hz range. The 
rms amplitudes are obtained by integrating the normal- 
ized PDS and corrected for b oth the background an d the 
ridge emission as described in lKalemci et al.l (|2006al) . All 
spectral and timing results are presented in Table [1] 

2.2. Swift Observations 

We also analyzed Swift X-ray Telescope (XRT) ob- 
servations conducted at the same time period with the 
RXTE observations. We found 12 observations carried 
out between MJD 55,622 and 55,647. We used them to- 



gether with the RXTE observations and looked for pres- 
ence of any spectral softening. 

We analyzed the XRT photon counting mode event 
data using XRTPIPELINE task provided in FTOOLS 
package. Pile up was an issue for the first four ob- 
servations whose count rate were greater than 1 c/s. 
To remove its effects, following the SWIFT SCIENCE 
DATA CENTER (SSDC) recommendations, we selected 
the source photons from a ring with an inner radius of 
5" and an outer radius of 40". For the rest of the ob- 
servations, the source photons were selected in a circular 
region with a radius of 40". The background photons 
were accumulated from a ring with an inner radius of 
70" and an outer radius of 100" centered at the source 
position. 

For the spectral analysis, only events with grades 
0-12 were selected. The auxiliary response files 
were created by XRTMKARF and corrected using 
the exposure maps, and the standard response matrix 
swxpc0tol2s6.20010101v013.rmf was used. We binned 
the energy spectra by fixing the number of counts per 
bin at 50. We fitted the spectra with a model that con- 
sists of photoabsorption and a power law in 0.6-8.0 keV 
band. In our initial spectral runs, we let the Njj free, and 
the resulting values were between (0.3 — 0.7) x 10 21 cm -2 . 
As we were not able to constrain the Nh, we performed a 
second run with Nh fixed at 0.5 x 10 22 cm -2 . The log of 
Swift observations and the spectral results are presented 
in Table □ 

2.3. SMARTS Observations 

The regular optical/infrare d observations wer e per- 
formed with the ANDICAM (jDePov et al.l I2003D cam- 
era on the SMARTS 1.3m telescope in V, I, J and H 
bands. The observations covered the outburst decay in 
daily basis between MJD 55,582 and 55,720. In this pa- 
per, we focus on the OIR light curves and evolution of 
the spectral energy distributions (SED). The deredden- 
ing of the observed magnitudes and their conversion to 
physical units were critical to create SEDs. For this pur- 
pose, we used the op tical extinction, Ay = 3.7 ± 0.3 
(jZdziarski et al .111 9981) together with the extinction laws 
given bv iCardelli et al.l (|1989h and lO'Donnelll (11994ft 
The same Ay was previously utilized inlCorbel fc Fender! 
I2002t ICoriat et al.lhooil IBuxton et~alll2012l for the SEP 
creation. For the details of the selection of Ay, dered- 
dening and the flux conversion procedures, we refer to 
IBuxton et aD (pOll ). 

3. RESULTS 

3.1. X-ray Evolution 

In Figure [TJ we present the evolution of the spectral 
and the temporal parameters. On MJP 55,594, dra- 
matic changes occurred in both the X-ray spectral and 
the temporal parameters. The rms amplitude of vari- 
ability jumped from 9.8% to 17.4% in one day. Both 
the disk flux and the temperature of the inner disk de- 
creased. The power law flux increased and the photon 
index started to harden. These changes in the evolu- 
tion of the parameters suggest a reshaping of the accre- 
tion dynamics and lead to a transition from the soft- 
intermediate state towards the hard state. In Figure [2j 
we show the hardness-intensity diagram and mark the 
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Observational Parameters Obtained From RXTE Data 
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14 


0.59 ± 


02 


4.14 


7.33 


<5.90 




01-02 


55,565.82 


0.53 


2 


46 


± 
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0.58 ± 
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8.71 
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<2.73 
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2 


27 
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0.56 ± 


01 


5.95 


4.37 
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03-00 
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01 
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0.56 ± 


01 
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55,580.61 


1.32 


2 


45 
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06 


0.55 ± 


02 


8.05 


3.62 


8.08 ± 0.53 
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55,582.70 


0.54 


2 


57 


± 





12 


0.54 ± 


03 


7.32 


3.09 


<7.50 
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55,585.94 


0.98 


2 


32 
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04 


0.56 ± 


01 


7.62 


3.20 


8.88 


± 0.19 


2.09 ± 0.06 


04-02 


55,586.49 


0.70 


2 


16 


± 





10 


0.55 ± 


03 


7.56 


2.96 


8.52 ± 0.73 




04-07 


55,587.50 


0.52 


2 


38 


± 





12 


0.54 ± 


03 


6.38 


2.35 


7.82 ± 1.00 




04-08 


55,588.55 


0.54 
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36 
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14 


0.53 ± 


03 


4.95 


2.53 


<10.94 




05-00 


55,589.20 


0.58 


2 


35 


± 





12 


0.52 ± 


03 


5.97 


2.30 


<9.90 
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55,590.43 


0.58 
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13 
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12 


0.51 ± 


04 


6.10 


1.87 


<10.61 
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55,591.61 


1.31 
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15 
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0.53 ± 


02 
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2.59 


8.63 ± 0.25 


1.77 ± 0.06 
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0.83 
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35 
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0.51 ± 


01 


6.25 


1.67 


11.30 ± 0.87 


1.72 ± 0.05 


05-02 


55,593.50 


1.55 
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li 
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08 


0.52 ± 


02 


5.89 


2.01 


9.83 ± 0.43 
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05 


7.31 


1.10 
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± 2.56 
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55,597.25 


0.98 
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01 


± 





03 


0.42 ± 


05 


9.27 


0.48 


19.38 


± 1.29 


1.03 ± 0.05 


06-01 


55,598.66 


1.09 


2 


11 


± 





06 


0.44 ± 


08 


7.94 


0.50 


25.18 


± 1.48 




06-02 


55,601.88 


1.80 


1 


89 


± 





03 


0.44 ± 


18 


8.35 


0.11 


22.94 


± 1.71 


1.39 ± 0.05 


07-00 


55,603.98 


0.74 


1 


77 


± 





04 






7.97 


0.00 


17.46 


± 1.62 




07-03 


55,604.89 


1.28 


1 


75 


± 





02 






7.61 


0.00 


26.16 


± 1.93 




07-01 


00,006.89 


1.05 


1 


67 


± 





04 






6.82 


0.00 


25.02 


± 1.12 




07-02 


55,607.76 


1.26 


1 


71 


± 





02 






6.05 


0.00 


23.04 


± 1.96 
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55,609.84 


1.24 


1 


66 
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02 






4.82 


0.00 
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0.00 


37.76 
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55,615.45 


0.85 


1 
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11-00 


55,632.05 
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77 


± 
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0.00 


44.79 


± 4.17 




11-02 


55,636.20 


0.59 


1 


68 


± 





10 






0.70 


0.00 


43.23 


± 5.35 




12-00 


55,638.73 


1.46 


1 


69 


± 





12 






0.59 


0.00 


41.03 


± 3.58 




12-01 


55,639.50 


1.52 


1 


86 


± 





14 






0.50 


0.00 


51.31 ± 10.72 
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1 


38 


± 
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0.31 


0.00 
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1 


72 


± 
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Observational Parameters Obtained From Swift Data 
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a Full observation ID is, 95409-01-Obs for the first observation, and 96409-01-Obs for the rest. 

b Modified Julian Date (JD-2,400,000.5) at the start of the observation. 

c Unabsorbed power law flux in the 3-25 keV band, in units of 10 — 10 ergs cm -2 s _1 . 

d Unabsorbed disk black-body in the 3-25 keV band, in units of 10 — 10 ergs cm -2 s — 1 . 

e Total rms amplitude of variability integrated over a range of O-oo Hz in the 3-30 keV band. 

f QPO centroid frequency 
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Figure 1. Evolution of (a) the total rms amplitude of variability 
in the 3-30 keV band, (b) the power law flux (filled circles) and 
the disk-blackbody flux (empty circles) in the 3-25 keV band in 
units of 10 — 10 ergs cm -2 s —1 , (c) J-band IR light curve, (d) the 
photon index T, (e) the inner disk temperature Tj n , (f) the ratio of 
the power law flux to the total flux in the 3-25 keV band. Dashed 
line indicates the time of state transition on MJD 55,594. Trian- 
gles and squares show the observations with type B and C QPOs, 
respectively. 



start of the state transition. Based on the extrapola- 
tion of the power law model, the X-ray luminosity at 
the start of the transition is Li_200 koV = 2.7 x 10 37 
ergs s" 1 or ~ 2% Ly.dd if we adopt a di stance of 8 kpc 
(iHvnes et all 120041 iZdziarski et aUl2004j) and a mass of 
10 M g (jHvnes et al.ll2003al iMunoz-Darias et aLll2008h . 



Before the transition, the energy spectra were soft with 
a mean photon index of 2.37. There was comparable con- 
tribution to the total X-ray flux from the disk and the 
power law components. The power law ratio (PLR, the 
ratio of the power law flux to the total flux in the 3-25 
keV band) was increasing from 0.40 to 0.75 level. This 
increase was due to the steady decrease in the disk flux. 
In addition to the decreasing disk flux, the inner disk 
temperature was also decreasing. The first detection of 
the X-ray variability occurred when the power law ratio 
reached 0.7. The rms amplitude of variability was less 
than 10%. Our an alysis confirms the type B QPO de- 
tections reported in lStiele et al.|[20Tll All these spectral 
properties indicate that the source was in transition from 
the so ft to the hard state, or simply in the intermediate 
state (jKalemci et al.ll200l . 

After the transition, the energy spectra became dom- 
inated by the power law component in six days. At the 
same time, the photon index hardened from 2.3 to 1.8. 
The power law flux increased, and remained at a higher 
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Figure 2. Start of the state transition and the approximate time 
of the OIR rcbrightcning arc marked in the hardness-intensity di- 
agram of the entire outburst. Observations analyzed in this paper 
are shown with filled circles (colored red in the online version). 
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Figure 3. Evolution of the photon index obtained from RXTE 
(circles) and Swift (diamonds) observations between MJD 55,611 
and 55,650. Best fit values at 1.53 and 1.70 constant levels for 
RXTE (dotted line) and Swift (solid line) only points, respectively. 

level than its soft-intermediate state level. Again during 
this six days the rms amplitude of variability increased 
to 25%. For two obse rvations, Type C Q POs (according 
to the classification in lMotta et al.l (|201ll )) were detected 
(see Figure Q] and Table [p. After MJD 55,605 the disk 
component was no more significant and no longer needed 
in the energy spectra. The power law flux started to de- 
cay and the photon index kept decreasing until it leveled 
off at 1.70. The rms amplitude of variability gradually 
increased to 50% in thirty days. 

3.2. No Evidence For Softening 

We also inspected the possible presence of softening of 
the spectra at low flux. In Figure [31 we plotted the evolu- 
tion of the photon index obtained from both the RXTE 
and Swift observations between MJD 55,611 and 55,650. 
The RXTE indices are systematically higher than the 
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Figure 4. RXTE/PCA X-ray, and SMARTS OIR light curves 
of GX 339-4. (a) X-ray flux in the 3-25 keV band in units of 
10 -10 ergs cm -2 s . Filled circles: power law flux (PL), empty 
circles: disk flux, (b-e) Undereddened H, J, I and V light curves in 
units of mjy. (Error bars are smaller than the plot symbols.) The 
dashed line indicates the X-ray state transition and the dotted line 
indicates the start time of the OIR rebrightcning. The solid lines 
show the baseline emission possibly originating in the disk. The 
arrows point the dates for which we constructed the SEDs (MJD 
55,609.84 and MJD 55,617.53). 

Swift indices, however both data sets are separately con- 
sistent with a flat evolution (T = 1.53 ±0.06, 1.70 ±0.01 
for Swift and RXTE, respectively). The reason for the 
systematic difference may be caused by the use of differ- 
ent energy bands in RXTE and Swift spectra. If the 
Galactic ridge emission is underestimated, the RXTE 
spectral indices would become slightly harder, but not 
enough to account for the entire difference. Regardless 
of the deviation between two data sets we conclude that 
the data suggests no evidence for the softening of the 
spectra between MJD 55,610 and 55,650. 

3.3. Light Curves 

In Figure IH we present the evolution of the power law 
and disk black body fluxes together with the OIR light 
curves obtained during the 2011 decay. The dashed line 
shows the time of state transition from soft-intermediate 
to hard-intermediate state whereas the dotted line shows 
the start time of the OIR rebrightening. To find the 
start of the rebrightening, we first formed a baseline 
that smoothly connects the fluxes before (between MJD 
55,590 - 55,604) and after (MJD 55,680 - 55,690) the flare 
as an exponential decay. We assumed that the physical 
origin of the rebrightening is separate from that of the 
baseline. We then fitted the rise of the rebrightening (be- 
tween MJD 55,608 - 55,612) in the infrared bands only 



with a straight line over the baseline. We used the in- 
frared, since the fluctuations are lower during the early 
part of the rebrightening compared to the optical bands. 
The start of the flare is defined as the date that the linear 
fit intersects zero, and it is MJD 55,607 ± 1 day. 

The evolution of the disk and the power law flux were 
described in § 13.11 As the source transitioned to the hard- 
intermediate state the power law flux increased, and the 
OIR rebrightening occurred when the spectrum was al- 
most its hardest (see also Figures [T] and [5]). There was 
a delay of ~12 days between the increase in the power 
law X-ray flux and the OIR rebrightcning. The rebright- 
ening started at a PCA flux of ^3^25 kcV=7-6 x 10~ 10 
ergs cm~ 2 s _1 , resulting in a bolomctric luminosity of 
£1-200 keV=l-8 x 10 37 ergs s _1 or - 1.4% L Edd . Note 
that such a delay has already been noticed for other black 
hole tr ansients, name ly XTE J1550 - 564 (Kalemci et al.l 
l2006bD . 4 U 1543-47 (IKalemci et all 120051) and also for 
GX 339-4 (|Coriat et al.ll2009D . 

The evolution in OIR in different bands are similar 
(see Figure |4]). A decay is followed by a ^70 days of 
rebrightcning that peaked around the same dates before 
reaching a constant level. The amount of brightening is, 
however, different among the bands. Ratio of the peak 
flux to the baseline flux decreases from H to V (max ^4.9 
to -3 on MJD 55,620). 

3.4. Evidence For Binary Period in The Optical Light 

Curves 

OIR light curves shown in Figure 0] fluctuate during 
the initial decay (between MJD 55,580 - 55,605) in a 
time scale of days. The fluctuations continue even on 
the rise and the peak of the rebrightcning. Moreover, 
some parts of the light curves seem to show regular mod- 
ulations. Therefore, we decided to search for per iodicity 
in all bands using the Lomb-Scargle algorithm ([Scarglel 
[1982] ). The initial decay (between MJD 55,580 - 55,605) 
did not provide a significant peak in the periodogram 
in any of the bands. Likewise, there is no evidence for 
periodicity in the light curves after MJD 55,630. 

On the other hand, we detected the known binary 
period of the system a t 1.77 days ()Hvnes et al.l l200"3ai 
iLevine fc Corbetl [2006) in the V, I and J band light 
curves between MJD 55,605 and 55,621, which is the rise 
and the peak of the rebr ightening. Using the method of 
I Home fc Baliunasl (|1986| ) the false alarm probabilities of 
the known period are estimated as 2.2xl0 -4 , 1.6xl0~ 4 , 
7.6x10-3 (3.69-cr, 3.78-cr, 2.78-cr) for the V, I and J 
bands, respectively. We performed a detailed fitting with 
the following function: 



2tt 

a + bt + ct 2 a + dtl + A sin(—t + 0) (1) 



where to is defined as the time from MJD 55,605, and F oir 
is in units of mjy. The polynomial part of the Eq. [1] rep- 
resents the continuum of the rebrightening whereas the 
periodic part represents the modulations with the binary 
period P = 1.77 days obtained from the periodogram. 
The model parameters obtained from the fits are given 
in Table [5] and the best model fits to the data are shown 
in Figure O The phases are consistent among all bands. 
Note that for the H band the periodic modulation is not 
necessary to fit the data, however the first few points are 
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consistent with the binary period if the modulation is 
included in the fit. 

3.5. SEDs 

We created spectral energy distributions to inspect the 
physical components and the mechanisms that produce 
the observed emission in both the intermediate state and 
the hard state. Below we will present our observational 
results by exploring the evolution of the OIR SED alone, 
and then together with the X-rays. 

3.5.1. OIR SEDs 

Figure [5] shows several SEDs of the total emission sam- 
pled from different stages of the OIR evolution. In order 
to quantitatively track the evolution, each SED was fit- 
ted with a power law. The evolution of the SED started 
with higher spectral slopes (all slopes are negative, and 
a represents the absolute value of the slopes throughout 
the paper) in the hard intermediate state (a = 1.33 on 
MJD 55,595) and continuously decreased (a = 0.76 on 
MJD 55,611) until the peak of the OIR rebrightening in 
the hard state. During the peak, the slope of the SED 
became flatter (a = 0.35 on MJD 55,618). As the OIR 
decayed, the slope started to increase again (a = 0.72 on 
MJD 55,650). At the end of the decay where the OIR 
emission decreased to a constant level in all bands, the 
spectral slope had a less steep value (a = 1.07 on MJD 
55,700) than it was in the hard intermediate state. The 
typical errors for the spectral slopes were around ±0.15 
at 1-cr, but note that the relative slopes are not affected 
by the dereddening process which dominates the error 
calculation. 

3.5.2. Broadband SEDs 

In order to understand the emission mechanism that 
gives rise to the rebrightening, we constructed two broad- 
band SEDs in Figure [7] with the total OIR (circles), the 
excess OIR (diamonds) , and the X-ray fluxes from MJD 
55,609.84 and 55,617.53 (rise and peak of rebrightening, 
see Figure 0]). The total OIR fluxes were calculated us- 
ing the polynomial part of Equation 1 (the solid curves 
in Figure [5]) This way we removed the effect of modu- 
lations with the binary period found in § 13.41 We then 
removed the baseline flux to find the SED of the excess. 
Wc calculated the errors on the flux densities by only con- 
sidering the calibration and interstellar extinction. The 
errors from the baseline and the rebrightening fits were 
ncgligablc. The excess OIR spectra yielded flat spectra 
with spectral slopes of 0.15 ± 0.15 for MJD 55,609.84 
and -0.05 ± 0.15 for MJD 55,617.53. We refer to the 
Discussion section for the explanation of the lines in Fig- 
ure [71 

4. DISCUSSION 

The changes observed in both the spectral and the tem- 
poral parameters imply a change in accretion dynam- 
ics. The most clear changes in both the X-ray spectral 
and temporal parameters occurred on MJD 55,594. We 
marked this date as the start of the transition from soft- 
intermediate state towards the hard state. Throughout 
the the decay phase of the outburst, GX 339-4 generally 
presented typical X-ray and OIR behaviour as observed 
from other black holes, and also of its previous outburst 
decays. 



4.1. Evolution In The X-ray Regime 

In the soft-intermediate state, the evolution in X-ray 
and OIR regimes was quite typical. Smoothly decay- 
ing disk temperature and disk flux may imply presence 
of a steadily cooling accretion disk. Variable power law 
flux with a variable soft photon index can be thought 
as arising from a variable corona which have some over- 
lap with the accretion disk that provides strong cooling. 
The emission mechanism that produces the power law 
com ponent is most probably the thermal Comptoniza- 
tion (jSunvaev fc Titarchuk|[l980T) . 

After the start of state transition, the disk temper- 
ature and the disk flux rapidly decayed and remained 
out of the PCA energy range in about ten days. As 
the disk is thought to be the source of seed photons 
in the Comptonization process, this is most often in- 
terpreted as the reduction in the number of soft seed 
photons. This results in a less cooling of the corona, 
and higher electron temperature which is consistent with 
the observed hardening of photon index. This simple 
Comptonization scenario could change or be modified 
with the onset of compact jet. Because the jet can ei- 
ther provide additional source of seed photons from non- 
thermal synchrotron emission at the base of jet (corona) 
and channel some portion of the accretion energy into its 
own power (jFender et al J 120031) . These processes could 
provide a sharp increase in the rms amplitude of vari- 
ability (|Hvnes et al.ll2003bl ; IRussell et alil2011bD a nd /or 



a softening of the X-ray spectrum (|Dincer et al.l 120081 : 
IRussell et al J [201(1 ). however no such sharp changes in 
the X-ray spectral and temporal properties were observed 
during rebrightening in the OIR. 

4.2. Origin Of OIR Emission On The Initial Decay 

At the initial decay (between MJD 55,585 and 55,605, 
sec Figure[3|, the OIR flux was most probably dominated 
by the thermal emission originating from the outer parts 
of the accretion disk as the secondary star is not expected 
to co ntribute to the OIR flux significantly (Shahbaz et al. 
120011 ). This is supported by the higher spectral slopes ob- 
served compared to the slopes during the rebrightening 
(see Figure [5]). Additionally, the spectral slope of 1.34 ± 
0.16 in this stage is lower than the expected spectral in- 
dex from the Rayleigh Jeans tail of the black body (F u oc 
v ) and s eems consistent with the X-ray irradiation of 
the d isk (jHvnesI 120051 iCoriat et al.l 120091: iBuxton et all 
I2012D . The strong variability in the light curve in the 
intermediate state also supports X-ray irradiation. Fur- 
thermore, recent optical and near-infrared spectroscopic 
observations at the soft state suggests that the irradia- 
tion might be enhanced due to a larger illumination area 
provided by a warped disk, and/or the winds launched 
at large radii of the accretion dis k that up-scatter pa rt 
of the X-ray and/or UV emission (jRahoui et all 12013) . 

4.3. Rebrightening Due To A Jet? 

The occurrence of OIR rebrightening in the hard 
state decay is similar to th ose observed in other 
black hole binaries 4 U 1543-47 (jBuxton fc Bailvn|[20p: 
Kalemci et all 120051) an d XTE J1550-56 4 dJain et~al 



20011: IRussell et al.ll2010D . XTE J1752-223 (|Russell et al 
2011bf ). Such a rebrightening has also been observed 



in 2003, 2005 and 2007 hard state decays of GX 339- 
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Table 2 

Model Parameters Derived From The Light Curves (for a fixed period P = 1.77 days) 



Parameters V I J H 

a (mjy) 0.229 ± 0.001 0.479 ± 0.007 0.918 ± 0.017 0.820 ± 0.050 

b (mjy d- 1 ) 0.047 ± 0.002 0.105 ± 0.004 0.200 ± 0.029 0.213 ± 0.026 

c (mjy d~ 2 ) -0.0029 ± 0.0003 -0.0031 ± 0.0006 0.017 ± 0.004 0.033 ± 0.004 

d (mJy) (7.23 ± 1.25) x 10~ 5 (-2.27 ± 2.62) x 10~ 5 -0.0011 ± 0.0002 -0.0021 ± 0.0001 

A (mjy) 0.061 ± 0.001 0.136 ± 0.003 0.216 ± 0.02 0.17 ± 0.02 

<j> (rad) 0.40 ± 0.02 0.48 ± 0.02 0.320 ± 0.079 0.232 ± 0.105 



4 (jCoriat et all 120091 : iBuxton et al.ll2012f ). In all these 
cases, the rebrigtening occurs at an X-ray luminosity 
range of 0.08 % to 2% Ls dd (Kalemci et al. in prep.). 4U 
1543- 47 (IPark et al.l2004D , XTE J1550-564 (ICorbel et al.l 
[2001 and XTE J1752-223 (|Miller- Jones et al.ll2011[ ) have 
shown radio revival sometime during the rebrightening, 
but the coverage for these cases were not adequate to 
describe radio behavior during the OIR rise. For 4U 
1543-47, XTE J1550-564 and XTE J1752-223, the SED 
of the excess during the rebrightening has indicated clear 
negative slope which is consistent with the optically thin 
synchrotron emission from a compact jet, and an extrap- 
olation from radio to OIR matches well with a flat or 
slightly inverted spectral slopes. GX 339-4 during the re- 
brightening in th e 2005 hard state decay showed similar 
radio-OIR SEDs (|Coriat et al.ll2009t) . Note that rebright- 
ening in the hard state decay of GX 339-4 in 2011 is also 
associated with radio revival (S. Corbel, private comm.), 
hence it is natural to assume that the rebrightening has 



jet origin. 

4.4. Understanding The Broadband SEDs 

We have produced two broadband SEDs from the rise 
and the peak of the rebrightening. Unlike 4U 1543-47 
and XTE J1550-564 for which a clear negative slope is 
present, our baseline subtracted SEDs indicate a flat to 
slightly inverted spectrum (see Figure [7]). Assuming that 
the excess is from the jet, we tried to obtain the break 
frequency for the change from an optically thick to an op- 
tically thin synchrotron emission, which is an important 
parameter in determining the base radius of the jet, ma g- 
netic field and minimum total jet power (|Fenderll2006l) . 

The flat SEDs may indicate that all OIR points lie 
on the optically thick part. In this case the break is 
around the V band, or at a higher frequency, which 
places a constraint on the slope of the optically thin syn- 
chrotron emission. A power law fit from V-band that 
passes through the softest X-ray band would give the 
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Figure 6. Dcrcddcncd OIR SEDs sampled from different stages 
of both intermediate and hard states. The numbers indicate the 
dates in the form of MJD-55,000. The dashed lines show the power 
law best fits to the data from individual days. 

spectral slopes of a = 0.72 ± 0.03 and a = 0.95 ± 0.03 
for MJD 55,609.84 and 55,617.53, respectively (the dot- 
dashed lines in Figure [7]) . These slopes are a lower limit 
to the possible absolute slopes that do not provide X- 
ray fluxes higher than the observed ones. The former is 
a typical spectral slope expected from an optically thin 
synchrotron jet emiss ion, but the latter is a bove the given 
range 0.6 < a < 0.8 ([Russell et aJJl20lTbh depending on 
the lepton energy distribution. If the spectral slop es as 
steep as 1.5 for XTE J1550-564 (IRussell et al J 12010) and 
XTE J1752-223 ([Russell et al.ll2011bD are possible at the 
beginning of the OIR rebrightening, then the particle dis- 
tributions inferred from SEDs would not be problematic. 

If the break is at the V band, however, the ratio 
of the minimum integrated synchrotron power to the 
X-ray luminosity is -Lj/Li_200 koV=0.18 and 0.95 for 
the SEDs on MJD 55,609.84 and MJD 55,617.53, re- 
spectively. For this calculation, we considered only 
the optically thick part of the SEDs and assummcd 
a conservative jet rad iative efficiency of 0.05 following 
ICorbel fc Fe"ndel[2002l The reported ratios for G X 339-4 
arc in between 0.05-0.10 (ICorbel fc Fender|[200l and the 
maximum valu e reported so far is for XTE J1118+480 
which is 0.20 ([Fender et al.l [200l . The ratio of 0.95 
is t oo high to be accoun ted by scale-invariant jet mod- 
els ([Markoff et all [20051) . The ratio deer eases down to 
0.32 for a break at the H band for this observation, and 
0.05 for the observation on MJD 55,609.84. The typi- 
cal break frequen cies reported so far f or GX 339-4 are in 
the mid-infrared (|Gandhi et al.l f201 If) to i nfrared range 
(ICorbel fc Fenderll2002t ICoriat et alj|2009t IHoman et all 
l2005HBuxton et al.H2012h~ 

On the other hand, if we place the break at H-band or 
at a lower frequency, we encounter a problem with a sim- 
ple post-shock emission model. In Figure [3 the dashed 
lines show the limiting cases of the power laws that can 



be fitted from H-band to softest X-ray band. The other 
three OIR bands remain significantly above the power 
law even when the errors are considered. This shows 
that the entire excess SED cannot be expla ined by opti- 
cally t hin synchrotron emission from a jet. iCoriat et al.1 
(|2009fj have studied SEDs from a similar stage of the 
hard state decay and claimed that significant reprocess- 
ing may be present during the OIR rebrightening. In 
this work, we show that even taking out the baseline 
emission that comes with reprocessing, there is still an 
excess. If a spectral break occurs at/below the H-band, 
it either mus t have two compon e nts (such as pre-sh ock 
synchrotron, iMarkoff et"aH 120031: IHoman et aTll2005[ ) or 
a simple flat to inverted spectrum breaking in the near 
or mid infrared is not a good description of the jet be- 
havior during its launch. Such an additional c ompo nent 
has recently been observed by iRahoui et al.l (j2012t ) . A 
full SED fitting (including radio fluxes) that takes into 
account the emission from different parts of the jet is re- 
quired to find the break frequency, and this is beyond 
the scope of this paper. 

We also extrapolated a power law (dotted lines) with 
spectral index obtained from X-ray spectral analysis 
down to the OIR band. In both days, the power law re- 
mained substantially below the OIR. This shows that the 
X-ray spectra cannot be explained by pure synchrotron 
emission without assuming a second emission component 
from the jet, at least during the early part of the OIR 
rebrightening. X-rays domi nated by dire c t jet emission 
scenario suc h as stated in iMaitra et al.l ([20091 ) is also 
dismissed in iBuxton et al.l (|2012f ). The svnchrotron-self 
Compton models are a lso able to explain the broadband 
natur e of jet emission ([Markoff et al.l [20051 ICoriat et al.l 
l2009f ). However a detailed analysis of such models are 
beyond the scope of this paper. 

4.5. On The Modulations Of the OIR Light Curves 

We have found that the OIR light curves on the rise of 
the rebrightening are modulated with the binary period 
of the system (~1.77 days). This can be explained with 
the X-ray irradiation of the secondary star by the hard 
photons originating from the corona. Then, we need to 
explain why we do not detect variations with the binary 
period before the rebrightening, and after the peak of the 
OIR light curve. Between MJD 55,580 and 55,605 (be- 
fore the rebrightening), the power- law flux is relatively 
high, however, the spectral indices are still high. During 
this time, the size of the corona may be relatively small, 
and therefore cannot illuminate the secondary star effec- 
tively while still effectively illuminating the warped disk. 
During this time, variations in the OIR flux from the 
disk is larger than the variations with the binary period. 
When the corona becomes larger, as indicated by the 
changes in the spectral index and the flux, it not only 
illuminates the secondary star effectively, but also allows 
the magnetic flux to travel i n close to the black hole, aid- 
ing in the launch of the jet (Meier 2001; Bcckwit h et al.l 
l2009f ). At the peak, the lack of modulations with the bi- 
nary period i s most likely due to the strongly variable jet 
synchrotron ([Rahoui et al.l 120121 ) dominating the emis- 
sion in all bands. After the peak, the X-ray power law 
flux decreases so much that it cannot produce signifi- 
cant irradiation on the surface of the secondary star. In 
Figure [4[ the decrease in the X-ray power law flux after 
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Figure 7. Broadband SEDs for MJD 55,609.84 and MJD 55,617.53. Black circles: total flux density of the source, diamonds: baseline 
subtracted flux density. Dotted line: power law fit extrapolations from X-rays, dashed line: power law fit from H-band to softest X-rays 
(3 keV), dashed-dotted lines: power law fit from V-band to softest X-rays (3 keV). The power law indices from the fits are shown in the 
legends. The errors on the slopes are calculated from the difference between the two slopes obtained by the fits from mid points of the OIR 
flux densities and the 1 cr above the mid points. 



MJD 55,607 is accompanied by a decrease in the ampli- 
tude of variations in V-band which is in favour of X-ray 
irradiation of secondary star. 
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